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ABSTRACT 

Using stellar models we investigate the relation between mass and the spectral type 
of the secondary star in low-mass short-period compact binaries such as cataclysmic 
variables and soft X-ray transients. Allowing for different mass transfer rates and dif- 
ferent system ages prior to mass transfer we find that the secondaries should populate 
a band in the spectral type vs. mass plane. The mass M ms of a ZAMS star with the 
same spectral type as the donor is effectively an upper limit to the donor mass. The 
lower mass limit for a given spectral type depends on the mixing length parameter. 
If this is large, there is no lower limit if the spectral type is later than K6. The band 
width decreases from O.4M0 at K6 to less than 0.2Mo at K0. If the mixing length 
parameter is small, there is no lower mass limit if the spectral type is later than M2, 
and the band width decreases from O.2M0 at M2 to less than 0.1M Q for types earlier 
than KO. 

We also point out an error in the method suggested by Beekman et al. (1997) to 
estimate the primary mass in a soft X-ray transient. 
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1 INTRODUCTION 

Reliable mass estimates for the stellar components of com- 
pact binaries, i.e. cataclysmic variables (CVs) and X-ray 
binaries (XBs), are notoriously difficult to obtain. Except 
for the few double-lined eclipsing systems a successful mass 
determination hinges on an estimate of the orbital inclina- 
tion and on subtleties in denning the effective dynamical 
centre of the emission from the donor star (see e.g. Smith & 
Dhillon 1998). 

In practice, masses in short-period CVs or XBs are of- 
ten estimated with the assumption that the donor star (the 
secondary) obeys a main-sequence mass-radius relation. A 
simple consequence of Roche geometry and Kepler's law is 
that in a semidetached binary the orbital period P and the 
mean density of the donor are related via 



r 2 — rm 2 , 
usually with a ■ 
Ph ~9m 2 . 



1, r ~ 1, so that (hi) reads 



(2) 



(3) 



7T12 



1/2 



(1) 



(e.g. King 1988). Here Ph denotes the period in hr, r 2 and 
m 2 the radius R 2 and mass M 2 of the donor star in solar 
units. The coefficient k 2 is a weak function of mass ratio 
q — M2/M1 (Mi is the primary mass), with k 2 ~ 8.86 for 
q Ja 1. If the secondary is in the core hydrogen burning phase 
its radius can be expressed as 



In short-period 'j-driven' compact binaries, where the 
orbital angular momentum losses drive mass transfer, the 
donor star can deviate from this standard mass-radius rela- 
tion mainly for the two following reasons, 
(i) The secondary experiences continuous mass loss which 
perturbs thermal equilibrium. The surface luminosity is no 
longer in balance with the luminosity generated in the core 
by nuclear reactions, and the difference causes the star to 
expand or contract. Mass-losing main-sequence stars are 
smaller than in thermal equilibrium when they are predom- 
inantly radiative, i.e. for m 2 <; 0.6, while they are oversized 
if they are predominantly convective, i.e. for m 2 £ 0.6 (e.g. 
Whyte & Eggleton 1983; Stehle et al. 1996; Baraffe & Kolb 
2000). 

(ii) The donor star could have evolved off the zero-age main- 
sequence (ZAMS) prior to mass transfer. Once mass transfer 
is under way the nuclear state at the onset of mass trans- 
fer is effectively frozen in. Once sufficiently low-mass, these 
donors mimic stars that are older than a Hubble time. Re- 
cently, Baraffe & Kolb (2000; see also Kolb & Baraffe 2000) 
argued that the observed late spectral type of CV secon- 
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daries provide strong evidence that a large fraction of CV 
secondaries have evolved off the ZAMS prior to mass trans- 
fer. Similarly, theoretical arguments suggest that the major- 
ity of short-period neutron-star XBs form with fairly mas- 
sive donors (m-z <; 1), allowing them to age prior to mass 
loss. The constant r in (0) can therefore vary significantly. 

Applying (H) to estimate the donor mass would usually 
give a value which is larger than the real mass. The donor is 
undermassive compared to a ZAMS star that fills the Roche 
lobe in the given binary. 

In the following we show that the effective temperature 
of the donor star — and hence its spectral type — is much 
less sensitive to the effects of mass loss and nuclear ageing 
than the radius. Using the latest generation of low-mass star 
models we explore theoretically to what extent the spectral 
type provides a reasonable estimate of the donor mass. In 
addition, we reconsider Beekman et al.'s (1997) attempt to 
place a lower limit on the black hole (or neutron star) mass 
in soft X-ray transients. 
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2 THE STELLAR CODE 

For our numerical experiments we employ the latest gener- 
ation of low-mass star and brown dwarf models by Baraffe 
et al. (1995, 1997, 1998, henceforth summarized as BCAH 
models) which represent a significant improvement in the 
quantitative description of stars with mass ^ 1 M© . 

The main strengths of the models are in two areas: 
the microphysics determining the equation of state in the 
stellar interior, and the non-grey atmosphere models which 
enter as the outer boundary condition. The equation of 
state (Saumon, Chabrier and Van Horn 1995) is specifically 
calculated for very low-mass stars, brown dwarfs and gi- 
ant planets. Recent much improved cool atmosphere models 
(Hauschildt et al. 1999; see also the review of Allard et al. 
1997) now provide realistic atmosphere profiles, which we 
use as the outer boundary condition, and synthetic spectra. 
Chabrier & Baraffe (1997) have shown that evolutionary 
models employing a grey atmosphere instead, e.g. the stan- 
dard Eddington approximation, overestimate the effective 
temperature for a given mass, and yield too large a mini- 
mum hydrogen burning mass. 

For a more detailed description of the code and the 
numerous tests against observations we refer to Baraffe & 
Kolb (2000), Kolb & Baraffe (1999; 2000), and references 
therein. 

We obtain the spectral type of a stellar model from its 
calculated colour (J — K) and the empirical SpT— (I — K) 
relation established by Beuermann et al. (1998). 

One of the main uncertainty in stellar models is the 
treatment of convection. Here we use the standard mixing 
length theory, with the ratio a = Z m i x / Hp of mixing length 
^mix and local pressure scale height Hp as a free parameter. 
Except for the specific case of the Sun, a cannot be cali- 
brated very well by a comparison between observations and 
models. The chief reason for this is that the age and chemical 
composition of observed stars is usually not known indepen- 
dently (see Baraffe et al. 1998; Baraffe 1999). According to 
such comparisons, the most likely range for a is between 1 
and 2. Recent 2-dim. hydrodynamical simulations (Ludwig, 
Freytag & Steffen 1999) also show a variation of a with effec- 



Figure 1. Spectral type versus secondary mass: dependence 
on the evolutionary state of the secondary at turn-on of mass 
transfer. Dots denote ZAMS models. Solid curve, close to the 
dots: donor unevolved (standard sequence). Solid curve, above 
the dots: donor moderately evolved (X c = 0.16; initial donor 
mass M2 t i = 1.2 Mq). Dashed: donor close to the terminal main 
sequence (X c = 4 X 10 — 4 , M^^i = 1.0 Mq). All sequences calcu- 
lated with low mass transfer rate (1.5 X 10 — 9 MQyr - 1 ). 



tive temperature T c h and surface gravity g between typically 
1 and 2. We therefore adopt this range of values to character- 
ize the present uncertainty due to convection. The spectral 
type of chemically inhomogeneous stars (m J; 0.6M©) de- 
pends on the choice of a, while unevolved, fully convective 
stellar models are insensitive to a. 



3 THE DONOR MASS-SPECTRAL TYPE 
RELATION 

The effective temperature T c h of low-mass main-sequence 
stars is rather insensitive to their degree of thermal disequi- 
librium. They behave just like giant stars on the Hayashi 
line, with nearly constant effective temperature. This is well 
understood in terms of the analytic Hayashi line theory for 
fully convective stars (e.g. Kippenhahn & Weigert 1990). 
Similarly, T e ff is not very sensitive to the nuclear age of the 
star. Single-star tracks in the HR diagram run almost ver- 
tical while on the lower main-sequence. 

Using BCAH models we quantify the deviation from 
equilibrium spectral type, as a measure of the surface tem- 
perature, for the donor stars in j-driven compact binaries 
with various evolutionary histories. We simplify the prob- 
lem by considering sequences with constant mass transfer 
rate M. To bracket the evolution of real systems where 
M presumably varies with orbital period we calculate se- 
quences with a high transfer rate and with a low transfer 
rate. Test sequences where we computed the transfer rate 
self-consistently according to simple orbital angular momen- 
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Figure 2. Spectral type versus secondary mass: dependence on 
the mass transfer rate. Dots denote ZAMS models. Thick lines 
evolved donor (X c = 4 X 10" 4 ); dashed: M = 1.5 X 10~ 9 M Q yr- 1 
solid: M = 5 X 10 — 8 MQyr -1 . Thin lines: unevolved donors 
dashed: M = 10~ 8 Moyr -1 ; solid: M = 10" 7 M yr -1 . 



Table 1. Spectral type and corresponding donor mass M2 along 
the sequences displayed in Fig. pi The mass in the second column 
(taken from the sequence with initial donor mass M2 % = 1.2 Mq, 
X c = 4 X 1CT 4 , M = 1.5 X lCr 9 M0yr-\ a = 2) defines an 
approximate lower limit to the mass of any donor with spectral 
type SpT, the mass in the third column (from the unevolved 
sequence with M 2 ,i = 1.0 Mq, M = lO _8 M0yr _1 , a = 1) an 
approximate upper limit. There is no lower limit for SpT later 
than K6. 



SpT 


M 2 /M 


M 2 /M 




lower limit 


upper limit 


G5 


0.83 


1.00 


G6 


0.79 


0.98 


G7 


0.75 


0.97 


G8 


0.72 


0.95 


G9 


0.68 


0.93 


KO 


0.65 


0.91 


Kl 


0.61 


0.88 


K2 


0.56 


0.85 


K3 


0.52 


0.83 


Kl 


0.47 


0.80 


K5 


0.42 


0.77 


K6 


0.35 


0.74 


K7 


- 


0.70 


MO 


_ 


0.66 


Ml 


- 


0.60 


AI2 


- 


0.50 


M3 


- 


0.38 


Ail 


- 


0.24 



turn loss descriptions give indeed the same qualitative and 
quantitative results. 

We first discuss sequences calculated with solar metal- 
licity and a — 1 (Figs, [ll and 0). 

In particular, in Fig. hi we show the standard sequence of 
Baraffe & Kolb (2000). Here mass transfer starts from an un- 
evolved (ZAMS) donor with mass 1 Mq, proceeds at a con- 
stant rate 1.5 x 10 -9 MQyr -1 , stops when the donor becomes 
fully convective (at mass 0.21 Mq), and resumes at the lower 
rate 5 x 10~ n Mgyr -1 once the donor has settled back into 
thermal equilibrium. This sequence fits the period gap in the 
framework of disrupted orbital braking (e.g. King 1988). In 
addition, in Fig. H we plot two other sequences that start 
from a 1 Mq ZAMS star (thin lines; henceforth "unevolved 
sequences"), one calculated with M — 10~ 8 Mgyr -1 , the 
other with 10 -7 Mgyr -1 . The mass transfer rates in j-driven 
compact binaries are usually estimated to be well ^ 10~ 8 
M yr _1 (see e.g. Warner 1995, his Fig. 9.8, for CVs; Chen 
et al. 1997 for XBs). A white dwarf system with M = 10~ 7 
Mgyr -1 would not appear as a CV at all, but presumably 
as a supersoft X-ray binary (e.g. van den Heuvel 1992; see 
also Kinget al. 2000). 

Fig. hi shows that the effect of thermal disequilibrium 
is negligible along the standard sequence; this hardly differs 
from that for systems with a ZAMS donor (i.e. in thermal 
equilibrium). If M is higher than in the standard sequence, 
the spectral type is slightly later for a given mass (Fig. 0). 
Even the extreme sequence with M — 10~ 7 Mgyr -1 never 
deviates by more than O.lMg for a given SpT from this 
relation (Fig. 0, thin solid line). In other words, unevolved 
donor stars follow an almost unique SpT-M2 relation, the 
one given by the ZAMS. 

This relation is much less well-defined when we include 
systems where the secondary has already burned a signifi- 
cant fraction of its hydrogen supply in the detached phase 
prior to mass transfer. Note that along these "evolved se- 
quences" the donor is still in the core hydrogen burning 
phase, otherwise the system would evolve off to longer peri- 
ods (see.g. Savonije & Pylyser 1988). 

Standard common envelope arguments predict that 
CVs form over the whole range of secondary masses, 0.1 S3 
m,2 Ss 1, hence the majority of CVs should have essentially 
unevolved donors (Politano 1996, de Kool 1992). Yet the 
generally late spectral type in systems with Ph ii 5 sug- 
gests that a fair fraction of CVs actually contain an evolved 
secondary (Baraffe & Kolb 2000; Kolb & Baraffe 2000; see 
also Ritter 2000) . Similarly, evolutionary considerations sug- 
gest that most short-period neutron-star XBs form with a 
somewhat evolved main-sequence donor with mass <; 1 Mq 
(King & Kolb 1997; Kalogera & Webbink 1998; Kalogera 
et al. 1998). This is mainly because systems with a more 
massive secondary stand a greater chance of surviving the 
supernova explosion and associated mass ejection at neutron 
star formation. 

It is therefore important to consider evolved sequences. 
Fig. Ill shows — in addition to the standard sequence dis- 
cussed above — a moderately evolved low M sequence (ini- 
tial central hydrogen content X c = 0.16, initial donor mass 
Mi,i = 1.2 M Q , M = 1.5 x 10~ 9 Mgyr" 1 ) and a significantly 
evolved, low M sequence (X c = 4 x 10" 4 , M 2 ,i = 1.2 M , 
M = 1.5 x 10 -9 Mgyr -1 ) where the donor was almost at 
the terminal main sequence at turn-on of mass transfer. 
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Figure 3. Expected populated domain in the spectral type versus 
secondary mass plane. The band denned by the unbroken curves 
correspond to a = 1.9, the band denned by the dashed curves to 
a = 1. (Thick linestyle: evolved sequences with X c = 4 X 10 — 4 , 
M 2> i = 1.2M@, M = 1.5xlO _9 M yr _1 ; thin linestyle: unevolved 
sequences with M = 10 — 8 M yr — - 1 -.) 



fined by Tab. w, even when the quoted errorbars are taken 
into account: U Gem, IP Peg, CN Ori and A0620-00. The 
donor masses of two of these systems have recently been re- 
vised to much smaller values. For U Gem (M4-M5) , Long & 
Gilliland (1999) give M 2 ~ 0.41M Q as preferred value. Given 
the uncertainty on the inclination we consider the remaining 
disagreement with the limits in Tab. nl as too marginal to 
deduce any conflict with theory. Similarly, Beekman et al. 
2000 find the new value M 2 ~ O.33±S;JgM for the donor 
mass of IP Peg. Taking into account the range of spectral 
type determinations (M2-M5; Ritter & Kolb 1998; Smith & 
Dhillon 1998) this system is now in agreement with Tab. hi 
The mass estimate for the discrepant system CN Ori is sus- 
pect as it is based on radial velocity curves exhibiting sig- 
nificant phase shifts (Smith & Dhillon 1998). Finally, the 
donor star in A0620-00 clearly is very undermassive and 
must be highly evolved. As the accretor is a black hole the 
initial donor mass could have been much larger than 1.5M© 
(the value considered here and representative for CVs and 
neutron-star LMXBs). 

If the observed sample is restricted to double-lined, 
eclipsing CVs where the secondary mass can be determined 
independent of the main-sequence assumption (pi) for the 
secondary, the limits in Tab. hi are acceptable in all but 
two cases (data from Ritter & Kolb, 1998, updated). The 
two cases where the agreement is only marginal are U Gem 
(see above) and EM Cyg, where the recently re-determined 
donor mass, 0.99 ± O.12M (North et al. 2000), appears to 
be somewhat too large for the spectral type K3. 



From Fig. H we see that the mass transfer rate has the 
same relative effect on evolved and unevolved sequences: the 
evolutionary track is shifted vertically downwards, to later 
spectral types. The tracks in heavy linestyle are evolved 
sequences, one with high mass transfer rate (solid line; 
M = 5 x 10 -8 Moyr" 1 ; X c = 4 x 10 -4 , M 2 ,i = 1.0 M Q ), the 
other with low mass transfer rate (dashed; M = 1.5 x 10 -9 
M yr"\ X c = 4xl0~ 4 , M 2 ,i = 1.2 M ; same as in Fig. []]). 

The highly evolved track with low transfer rate (M = 
1.5 x lO _9 M yr _1 ) and the unevolved track with high trans- 
fer rate (M = lO _8 M yr _1 ) should bracket the location of 
compact binary secondaries. 

Fig. ra shows the effect of the mixing length parameter 
q on these limiting curves in the SpT vs. mass diagram. Un- 
evolved tracks become undistinguishable for predominantly 
convective donors (M 2 ^ O.6M ), while the most evolved 
sequences are quite sensitive to a. Table hi gives lower and 
upper limits to M 2 for a given spectral type. The range re- 
flected in this table is to some extent due to our ignorance 
of the appropriate mixing length parameter. 

The metallicity Z of the donor can affect the mass- 
spectral type relationship as well. Lower Z gives a lower 
mass for a given spectral type. However, as stressed by 
Beuermann et al. (1998), there is no evidence for low metal- 
licity secondaries in the currently observed sample of CVs. 
Hence we restrict our analysis to systems with solar metal- 
licity (Z = 0.02). 

The observations compiled by Smith & Dhillon (1998), 
when taken at face value, are broadly consistent with the 
band shown in Fig. H (see their Fig. 3). They do not, how- 
ever, help to constrain the mixing length parameter. Four 
out of 22 systems in their sample fall outside the region de- 



4 LIMITS ON THE PRIMARY MASS 

This work originally started as an attempt to improve on a 
a paper by Beekman et al. (1997: hereafter B97) who con- 
sidered limits to the primary mass in soft X-ray transients. 
Although it turned out that our results do not significantly 
improve these limits, we wish to comment on the method 
suggested by B97 as a whole, and point out an error in their 
analysis, which so far has not been corrected in the litera- 
ture. 

B97 made use of an inequality derived by King, Kolb 
& Burderi 1996 (hereafter KKB), which expresses the con- 
dition for a short-period low-mass X-ray binary to have an 
unstable accretion disc and thus appear as a transient. Un- 
fortunately B97 misunderstood the definition of the quantity 
77i2 in this inequality. 7712 is defined by KKB as the ratio of 
the secondary mass 777,2 to that of the ZAMS star, m lna (p) 
(in solar units), which would fill its Roche lobe in a binary 
of the same orbital period, 



7772 



7712 



*o>y 



(4) 



As is obvious from (|l|) the relevant main-sequence star has 
the same mean density as the secondary; hence the notation 
mms(p)- B97 confused rh 2 with the ratio m 2 /r7i ms (7?), where 
777ms (R) is the mass (in solar units) of a ZAMS star that 
has the same absolute radius R2 as the secondary. It can be 
shown that 



.(«) 



7712 
I 2/3 ' 



(•») 
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so that a lower bound on mz (which is required to obtain a 
primary mass limit) is 



7712 > 



(>■! 



-i 3/2 



Pms (Tgf 



(6) 



Here pi is the mean density of the secondary, (given di- 
rectly by the binary period, see Eq. (ph, and p m s(T c ff) is 
the mean density of the main-sequence star with the same 
spectral type as the secondary. This is weaker by a fac- 
tor [p2/pms(T e ff)] compared to B97's limit, which in turn 
weakens the limits on Mi. 

In addition to this obvious mistake, we caution that the 
limitations and uncertainties inherent in using the disc insta- 
bility criterion are probably too severe to give a reliable mass 
limit. This criterion assumes that the binary transfers mass 
as a result of angular momentum losses ('j-driven') rather 
than nuclear expansion of the secondary ('n-driven'), but 
allows for the fact that this star may be somewhat nuclear- 
evolved before mass transfer begins. It is therefore only ap- 
propriate for short orbital periods. Although the disc insta- 
bility criterion correctly represents the statistical properties 
of the observed sample of soft X-ray transients, its applica- 
tion to individual systems is suspect for at least two reasons. 

1. The criterion is based on a comparison of the secular 
mean mass transfer rate driven by magnetic braking and the 
critical mass transfer rate for disc stability. Hence uncertain- 
ties in the magnetic braking strength, in particular at long 
orbital periods, and in the critical rate (see KKB; King, Kolb 
& Szuszkiewicz 1997; Dubus et al. 1999) propagate directly 
into the mass limit. 

2. It is not clear if a given system is j-driven, in partic- 
ular if the period is ^ 10 h. The ultimate fate of a system as 
j- or n-driven is determined by its orbital period Pi when 
mass transfer begins. If this exceeds a certain 'bifurcation 
period' Pbn (Pylyser & Savonije, 1988, 1989) the binary 
eventually becomes a diverging n-driven system, while if 
Pi < Pbit it ends up as j-driven. Uncertainties in for ex- 
ample the strength of magnetic braking at periods <; 12 hr 
make Pbif correspondingly uncertain. Further, even given a 
definite value of Pbif, n-driven systems born near this pe- 
riod initially evolve to periods below it (Pylyser & Savonije, 
1988, 1989) before nuclear evolution eventually wins out and 
drives the system to long periods. This comment is partic- 
ularly relevant for the neutron-star transients Cen X-4 and 
Aql X-l. 



5 SUMMARY 

We have applied up-to-date stellar models of low-mass 
main-sequence stars to investigate the mass-spectral type 
relation if these stars experience mass loss and are some- 
what nuclear-evolved. 

Unevolved stars follow an almost unique SpT-A/2 rela- 
tion, whatever the mass transfer rate. Nuclear-evolved stars 
have a somewhat earlier spectral type than a ZAMS star 
with the same mass. The spectral type is earlier if the star 
is more evolved. 

Accurate mass predictions for a given spectral type are 
hampered by the present uncertainty of the mixing length 
parameter. The mass M ms of a ZAMS star with the same 
spectral type as the donor is effectively an upper limit to the 



donor mass. The lower mass limit for a given spectral type 
depends on the mixing length parameter. If this is large, 
there is no lower limit if the spectral type is later than K6. 
The band width decreases from O.4M0 at K6 to less than 
0.2Mq at K0. If the mixing length parameter is small, there 
is no lower mass limit if the spectral type is later than M2, 
and the band width decreases from O.2M0 at M2 to less 
than 0.1M© for types earlier than K0. 

A logical next step is to consider the evolution of ther- 
mally unstable systems that start mass transfer from a fairly 
massive donor (M2 ^ 2Mq ) and reappear as j-driven short- 
period systems after reversal of the mass ratio. Detailed 
calculations of such systems are underway (Schenker et al. 
2000, in preparation). These will allow us to assess how dif- 
ferent descendants of thermal-timescale mass transfer are 
from the systems considered here, which had a more mod- 
erate initial mass (< 1.5M©) and were j-driven from the 
start. 

The results presented here are rather robust and in- 
evitable predictions of standard stellar structure and evo- 
lution theory that has been successfully tested against ob- 
servations of isolated low-mass stars. It is therefore highly 
desirable to test them by observational means, i.e. by a re- 
liable donor mass determination in systems with spectral 
type in the indicated range. 
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